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We present an embedding scheme for periodic systems that facilitates the treatment of the phys¬ 
ically important part (here the unit cell) with advanced electronic-structure methods, that are 
computationally too expensive for periodic systems. The rest of the periodic system is treated with 
computationally less demanding approaches, e.g., Kohn-Sham density-functional theory, in a self- 
consistent manner. Our scheme is based on the concept of dynamical mean-field theory (DMFT) 
formulated in terms of Green’s functions. In contrast to the original DMFT formulation for corre¬ 
lated model Hamiltonians, we here consider the unit cell as local embedded cluster in a first-principles 
way, that includes all electronic degrees of freedom. Our real-space dynamical mean-field embedding 
(RDMFE) scheme features two nested Dyson equations, one for the embedded cluster and another 
for the periodic surrounding. The total energy is computed from the resulting Green’s functions. 

The performance of our scheme is demonstrated by treating the embedded region with hybrid func¬ 
tionals and many-body perturbation theory in the GW approach for simple bulk systems. The total 
energy and the density of states converge rapidly with respect to the computational parameters and 
approach their bulk limit with increasing cluster (i.e., unit cell) size. 

PACS numbers: Valid PACS appear here 


I. INTRODUCTION 

Density-functional theory (DFT) has become a widely 
applied electronic-structure theory method due to the 
balance between accuracy and computational efficiency 
of local and semi local approximations such as the 
local-density (LDA) and generalized gradient approxi¬ 
mations (GGA). However, LDA and GGAs suffer from 
certain intrinsic limitations such as the self-interaction 
errorAi^, the absence of the derivative discontinuity in 
the exchange-correlation potentialAi^, the lack of long- 
range van der Waals interactions^^— , and the absence of 
image effectsiSrJ^. These shortcomings limit the predic¬ 
tive power of LDA and GGAs, in particular for localised 
electrons as found in d- or /-electron systemsi^^— or for 
adsorbates and surfaces— Furthermore, DFT is in¬ 
herently a ground-state method and therefore of limited 
applicability for excited states and spectra. More ad¬ 
vanced electronic-structure methods that overcome one 
or several of the mentioned shortcomings exist, but they 
are typically computationally much more demanding and 
thus limited to small systems sizes or a subset of elec¬ 
tronic degrees of freedom. To overcome the efficiency- 
accuracy conundrum, much effort has been devoted to 
combine the best of both worlds, that is to merge local 
and semilocal DFT approximations (DFA) with advanced 
electronic methods^^^— . 

We here advocate the concept of embedding. In this 
divide and conquer approach, the full system is divided 
into two parts: a small embedded region, which is treated 
with advanced, computationally demanding approaches, 
and an embedding environment that is treated with com¬ 
putationally more efficient approaches. A schematic il¬ 


lustration of the embedding concept is shown in Fig. [T] 



FIG. 1. Schematic depiction of the embedding concept. The 
embedded region (red) is treated with a more accurate theory, 
whereas the surrounding (grey) is calculated with less accu¬ 
rate and thus computationally less expensive theories. The 
key challenge of our and all previous embedding approaches 
is the appropriate treatment of the red/grey boundary. 

Following this general principle, various embedding 
schemes have been developed in the past^“— . They 
differ in scope (i.e., area of application), on how they 
treat the coupling between the embedded region and 
the surrounding, and in the approaches used to de¬ 
scribe the two regions. In (bio)chemistry, for example, 
one of the most popular embedding schemes combines 
quantum mechanics (QM) and classical molecular me¬ 
chanics (MM). The embedded region is treated quantum 
mechanically and the surrounding by MMi^“— In sur¬ 
face science, fully quantum mechanical schemes are more 
prevalent, e.g., for the description of surface adsorbates. 
They divide space into regions for advanced and less ad- 
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vanced electronic-structure approaches and differ mostly 
on how these two regions are coupled, e.g. through max¬ 
imal exchange overlap^^, density embedding^S or clus¬ 
ter extrapolation^SiSi. In solid-state physics, dynam¬ 
ical mean-field theory22ri2i offers a natural embedding 
framework by mapping an infinite, correlated lattice 
model into an impurity model^ immersed into a self- 
consistently determined mean-field bath. When DFA 
is chosen as the mean field, DMFT becomes material 
specifio^i^i^. The embedding is achieved by means 
of Green’s functions facilitating the calculation of spec¬ 
tra, band structures, but also phase diagrams. Recently, 
Zgid and Chait^ proposed to use DMFT as an embed¬ 
ding framework for quantum-chemistry approaches such 
as the configuration-interaction (Cl) method. They since 
proposed a simplified DMFT scheme based on density- 
matrix embedding to access static properties (e.g., the 
ground-state energy and its derivatives)^. However, at 
present all DMFT approaches use a down-folding proce¬ 
dure to a low energy subspace, which is treated on the 
model-Hamiltonian level. 

We here extend the DMFT concept to couple two first- 
principles regions. Our idea is similar to that of Zgid and 
Chan^^, but we explore the possibility of using DMFT as 
a general embedding scheme for advanced first-principles 
electronic-structure methods. These can be advanced 
DFT exchange-correlation functionals or excited-state 
methods based on the GW approach^, which are still 
computationally too expensive for large-scale systems. 
The difference to previous DMFT schemes is that we 
treat the unit cell as the local, embedded region, that is 
coupled to the rest of the periodic system via the DMFT 
framework. The advantage of our approach is that all 
electrons in the embedded region are treated on the same 
quantum mechanical level, which removes the arbitrari¬ 
ness in the definition of the down-folded subspace and 
does not require any double counting corrections. Fur¬ 
thermore, our approach includes non-local interactions in 
the embedded region, whose size can be systematically 
converged to the thermodynamic limit. 

We here present the concept of our real-space dynam¬ 
ical mean-field embedding (RDMFE) approach and its 
implementation in the all-electron Fritz Haber Institute 
ab initio molecular simulations (FHI-aims) code^“— . 
We first benchmark it for hybrid density functionals 
for which we have a periodic reference^. Then we 
apply our scheme to the GW approach. Due to its 
intrinsic self-consistency, our RDMFE approach yields 
a self-consistent GW solution, which is a much cov¬ 
eted approach for solids right now. While fully self- 
consistent GW implementations for molecules are slowly 
emerging^^— , we are only aware of one recent implemen¬ 
tation of fully self-consistent GW for solids^i^, which 
is, however, limited to small unit cells due to its com¬ 
putational expense. An alternative, approximate way 
to achieve self-consistency within GW is the so-called 
quasiparticle self-consistent GW (QPscGW) scheme^i^ 
which was recently been widely applied to solids. We here 


present self-consistent GW spectra and ground-state en¬ 
ergies for simple solids. 

The rest of the paper is organized as follows: The de¬ 
tailed formalism of our Green’s-function based embed¬ 
ding scheme is derived in Section |IT1 Section HVl presents 
benchmark results including both total energies and band 
structures for simple bulk systems. In Section 0 we con¬ 
trast our approach with the aforementioned embedding 
schemes. Section EH concludes the paper. 


II. SELF-CONSISTENT GREEN’S FUNCTION 
EMBEDDING IN REAL SPACE 

A. General concept 

In its original formulation^^ DMFT is a Green’s func¬ 
tion method for correlated model Hamiltonians that uses 
the locality of the electronic interaction to embed a local 
on-site region of the Hubbard lattice - typically a single 
d- or /-level - into a periodic electronic bath defining 
a self-consistent scheme. Treating the on-site region lo¬ 
cally facilitates the use of computationally very demand¬ 
ing and at the same time very accurate methods such 
as continuous time quantum Monte-Caro^, direct diag- 
onalization or renormalization group techniques^. The 
localized region is then coupled through a hybridization 
self-energy to the surrounding electronic bath, which is 
treated with computationally more efficient methods. In 
the last few years, DMFT has proven to be very success¬ 
ful in describing the spectral properties of solids with 
localised d- and /-states^^—. However, present versions 
of DMFT also suffer from certain shortcomings. For ex¬ 
ample, a double counting problem arises when local or 
semilocal DFA is used to calculate the electronic sur¬ 
rounding, because the DFA contribution of the localised 
manifold, which would have to be subtracted in order 
to not count it twice, is not known^. Different dou¬ 
ble counting schemes can give vastly different answers 
in DFA-I-DMFT— i^. Last, but not least, the choice of 
the local manifold can have a considerable impact on the 
result, in particular when the localised states hybridise 
with delocalised states in the system^i^. In principle, 
one could use renormalisation groups to derive the most 
suitable low energy Hamiltonian for the localised region 
systematically. However, in practice this is not tractable 
computationally. 

In this work, we use the concept of DMFT and for¬ 
mulate it as a Green’s function embedding scheme for 
the first-principles Hamiltonian. In contrast to the orig¬ 
inal DMFT formulation for correlated model Hamilto¬ 
nians, we consider the unit cell (or any computational 
supercell that can span the whole space when periodi¬ 
cally repeated) as embedded cluster in an first-principles 
way (see Fig. [2]). Our approach therefore includes all 
electronic degrees of freedom and does not require any 
downfolding. Our approach also permits the charge flow 
from one region to the other and therefore naturally in- 
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corporates the boundary between the two regions. No 
special treatment is necessary for atoms on the bound¬ 
ary, nor is it a problem with the boundary cuts covalent 
bonds. 


B. Embedding scheme based on DMFT 

1. Green’s function in a non-orthogonal basis set 


2. The “on-site” Green’s function for a periodic system 

The k-dependent Green’s function and self-energy in 
Eq. (|T|) can be Fourier-transformed to real-space, 

l.BZ 

Gif (R, - ^ e*(^--^^)''^Glf (k, 

k 

l.BZ 

Ey(R, - = - -V (k,ia;) 

i-^l.BZ , 
k 

(5) 


The embedding framework of DMFT is most conve¬ 
niently formulated in terms of Green’s functions. In a 
finite (and generally nonorthogonal) basis set the 

Green’s function G(r, r',zu;) can be expanded as, 

G{r,r',iu}) , (1) 

i,3 


where Gy {iuj) is the matrix form of the Green’s func¬ 
tion. Here we use the Green’s function on the imaginary 
frequency axis for computational convenience and with¬ 
out loss of generality. For a non-interacting Hamiltonian 
Rf = (e.g., the Kohn-Sham Hamiltonian), 

the corresponding non-interacting Green’s function G'^ 
in its matrix form satisfies 




where Sij = {(j)i\4)j) is the overlap matrix of the ba¬ 
sis functions, and n is the chemical potential. Using 
the Dyson equation that connects the non-interacting 
Green’s function G° with the fully interacting one (G), 
we obtain 

-I- yL)Sik - H'^k - ^ikiiu}))Gkj{iui) = Sij , (3) 

k 

where is the electronic self-energy. For periodic 

systems, the Hamiltonian and the Green’s functions are 
characterized by a Bloch k-vector in the first Brillouin 
zone of reciprocal space. Equation [3] thus becomes 


^((*o.+/r)5.fe(k)-R'i(k)-E,fe(k,*cc))Gif (k,*..) = % , 

k 

(4) 

with the lattice Green function G^*’'*(k, zw). Our imple¬ 
mentation is based on the all-electron FHI-aims code 
package?^ which uses numerical atom-centered orbitals 
(NAOs) as its basic functions. The basis functions {cfi} 
will thus be NAOs in our work. 


where Ri and Rj are Bravais lattice vectors denoting the 
unit cells in which the basis functions i and j are located. 
A^ibz is the number of k-points in the first Brillouin zone 
(l.BZ). The concept of DMFT is based on the fact that 
the lattice self-energy becomes local, or k-independent, in 
infinite dimension {D = oo)i^ For a crystal with trans¬ 
lational symmetry this implies 

Sy (R, - R„ zee) = El°^(iw),5R..R, . (6) 

Thus the self-energy is non-zero only if the two basis func¬ 
tions originate from the same unit cell. We call this the 
local (loc) or “on-site” self-energy, following the terminol¬ 
ogy of the model-Hamiltonian studies. In this limit, the 
whole periodic system can be mapped onto an effective 
impurity model of a local unit cell dynamically coupled 
to an effective “external” potential arising from the rest 
of the crystal. 

The first step to establish this mapping is to define 
the “on-site” Green’s function, i.e., Gy (Ri —Rj, zw) with 
Ri=Rj. Using the locality of the self-energy and Eq. o, 
we obtain the following expression for the onsite Green’s 
function. 


l.BZ 


g: 


on-site 


(iuj) = 


W 


l.BZ 


^Glf(k,z..) = 


l.BZ 


N- 


l.BZ 


^ [{IU3 + - R°(k) - Ei-(zc^)] . (7) 


In the DMFT context this equation is also known as the 
k-integrated Dyson equation. So far we have not speci¬ 
fied . In our embedding scheme, the environment is 
treated by KS-DFA in the LDA or the Perdew-Burke- 
Ernzerhof (PBE)^ GGA. A natural choice of is thus 
the KS-Hamiltonian i7^®(k) within LDA or GGA, that 
contains the kinetic-energy operator, the external poten¬ 
tial (zzext), the Hartree potential (ur), and the exchange- 
correlation (XG) potential (uxc) 

H^^{k) = -b z;ext(k) -b z;H(k) -b z;xc(k). (8) 

Next, we need to define E*°'’(za;) in Eq. ([7]). If we 
start from iJ^®(k), the “on-site” self-energy becomes 
the difference between the dynamic, complex many-body 
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FIG. 2. The DMFT embedding concept for a Si unit cell. The atoms in the unit cell (red region) constitute the embedded 
sub-manifold. Each unit cell of the periodic system is treated as a localized region, i.e. only local interactions are treated. 
The unit cells are coupled to the rest of the system via the hybridization self-energy (green arrow). 


exchange-correlation self-energy Exc(k, tw) and the KS 
XC potential, i.e., 

l.BZ 

E'-(*..) = --^ [Sxc(k,*t.) - t;^g(k)] 

-/Vi.BZ , 

k 

= ^xc(*^) ~ "^xc ■ (9) 

Using Eqs. 0 and ([9]), we finally obtain 
G°"-“‘®(iw) = 

l.BZ 

-- ^ [{^0J + ^,)S{k) - ijKS(k) _ S^oc ^ 

( 10 ) 

Our scheme is thus free from any double-counting ambi¬ 
guities, because the DFA XC-contribution that has to be 
subtracted is uniquely defined. 


3. Embedded Green’s function 

In the DMFT formalism, a periodic system is viewed 
as a periodically repeated cluster (here the unit or super¬ 
cell) dynamically embedded into a self-consistently deter¬ 
mined environment. The coupling between the embedded 
subsystem and its surrounding environment is described 
by a so-called bath Green’s function Q{iuj), connecting 
the Green’s function of the embedded cluster 
and the local self-energy via 

. ( 11 ) 

Here the local self-energy E^°'’(iaj) is the same as in¬ 

troduced in Fq. (0. The self-consistency condition of 
DMFT requires that the Green’s function of the embed¬ 
ded cluster G®™'’’(*a;) equals the on-site Green’s function 
as given in Fq. 0, 

G®”‘’(ia;) = G™-"‘*®(*w). (12) 

Alternatively, one can also use a so-called hybridiza¬ 
tion function A{iuj) to describe the coupling between the 


embedded cluster and its environment, which provides a 
more intuitive picture. A(zu;) is closely related to the 
bath Green’s function (/(fw), 

= (zw + ^l)S - - A(zw). (13) 

In Fq. (TT^ ^cluster jg Hamiltonian of the bare cluster 
describing the non-interacting unit cell i.e., without the 
"^xc contribution and without the presence of the other 
atoms from neighboring unit cells (see Fig. HD. This cor¬ 
responds to the “on-site” term of the Hamiltonian of the 
periodic system, and in practice can be conveniently ob¬ 
tained from the k-dependent Hamiltonian, 

.. l.BZ 

^cluster ^ ^ [^KS(k) _ ^KS(k)] . (M) 

tVi.BZ , 
k 

Using eqs. (HID-dnD, we obtain the following expression 
for the Green’s function of the embedded cluster 




j ij 


{iio + ^i)S - - S'°“[G®“‘’](ia;) - A(za;) 


j ij 


(15) 


Here we have explicitly indicated that the local self¬ 
energy is a functional of the embedded Green’s function. 
Thus eq. m has to be solved self-consistently, which 
corresponds to the inner loop of Fig. 2] The functional 
dependence of E'°'’(ia;) on G®“*’(iw) is given by the ac¬ 
tual approximation for the localized region, which will be 
the topic of next section. However, already here we see 
that our RDMFE approach lends itself to those advanced 
electronic-structure methods that can be expressed by 
(self-consistent) Green’s functions. 

Another point we would like to emphasize is the choice 
of the cluster overlap matrix S in eq. m- We found 
when updating the chemical potential of the cluster in 
the inner loop that we needed to define the cluster over- 
r l.BZ 1 


lap matrix as S 


1 

Ni.BZ 


E s-Hk) 


(and not sim- 
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l.BZ 

ply S = 'S'(k)) as done by Kotliar et alM to 

k 

enforce the correct asymptotic behavior of A(iw), i.e. 
lim A(a;)^ 0 . 


C. The local self-energy 

So far we had not specified the approximation for 
the local self-energy in eq. dTSl) . In our scheme the 
choice for yjioc quite flexible. In other 

words, we could use any approximation that goes be¬ 
yond LDA and GGAs. However, our framework lends it¬ 
self to Green’s-function-based approaches. This includes 
density-matrix and density-based approaches, because 
both quantities can easily be extracted from the Green’s 
function. Below we report on two different examples, 
namely hybrid density functionals that mix a fraction 
of exact-exchange with GGA semi-local exchange^^— 
and the GW approximation.— In practice, we could also 
go beyond GW, e.g., by including the screened second- 
order exchange (SOSEX) self-energy that was developed 
recently— 

We here use the PBE hybrid functional family 
(PBEh)^, whose most prominent functional is PBEO— 
We will also use the short-ranged range-separated hy¬ 
brid functional family by Heyd, Scuseria and Ernzerhof 
(HSE)— In PBEh the local self-energy in eq. (0 is given 
by 


^PBEh(Q^) — + (1 ~ + ^c"^] ~ "^XC 

= a . (16) 

In Eq. (ITSl) . is the “on-site” part of the GGA ex¬ 
change, and is the exact-exchange matrix given by 

= ( 17 ) 

k,l 


with SR and LR denoting the short and long-range part, 
respectively. If we now replace Ex'^’^^( 7 ) by and 

introduce the a parameter again, the local HSE self¬ 
energy assumes the following form 

E1?§e(«, 7) = a (Sf’‘°^(7) - (20) 

Furthermore, we employ the GW approximation for 
the local self-energy. Here, the computation of the GW 
self-energy for a given input embedded Green’s function 
follows the self-consistent GW implementation for finite 
systems in FHI-aims.— On the imaginary time axis, 
the GW self-energy for the embedded cluster is obtained 
as 


( 21 ) 

Here v indices refer to the auxiliary basis set used to 
expand the screened Goulomb interaction W^°^ in the RI 
approac h"^^d^ . Furthermore are the 3-index coeffi¬ 
cients obtained as, 




( 22 ) 


where 


and 


{ik\fi) = j drdr' 


,0i(r)((ife(r)P^(r') 


= / drdr' 


,P^ir)P,^ 


(23) 


(24) 


with {P^(r)} being the auxiliary basis functions. For W 
we thus obtain 

= E - n'-(zcc)]((i (25) 


where {ik\lj) are two-electron four-orbital integrals, and 
is the density matrix of the embedded cluster which 
can be obtained from the embedded Green’s function 

= (18) 

The two-electron Goulomb repulsion integrals are eval¬ 
uated using the resolution of identity (RI) technique in 
FHI-aims as documented in Ref. 113 The PBEO func¬ 
tional is obtained for a=0.25.— 

The extension to an HSE type self-energy is straight¬ 
forward. In HSE, a range-separation parameter is intro¬ 
duced that cuts off the exact-exchange contribution at 
long distances. The range is controlled via the screen¬ 
ing parameter 7 so that the local exchange self-energy 
becomes 

E^^( 7 ) = E^^'®’"( 7 ) + E^“’^^( 7 ), (19) 


where n'°‘’(za;) the irreducible polarisability, whose 
Fourier transform in the time domain is directly deter¬ 
mined by the embedded Green’s function 

n;r;(*r) = -z E (26) 

ijlm 

D. The self-consistency loops 

In our formalism eqs. (US and dig or m define 
an additional inner self-consistency loop for the local 
self-energy as depicted in Fig. H) Good convergence is 
achieved by a linear mixing 

E)(>;i = AE)(>“ + (l-A)E)(>fi, (27) 

with a mixing parameter A = 0.5. More advanced mix¬ 
ing schemes could be implemented as well, but we found 
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that linear mixing works well for the examples presented 
in this work. When the inner-loop reaches convergence 
we feed the resulting back into the on-site GF and 
iterate the main-loop further using the same mixing as 
for the inner-loop. 

Finally it is worth mentioning, that the on-site Green’s 
function as defined in eq. dT]) requires that our 1]'°'’ in the 
on-site Green’s function in the 0-th iteration should be 

Figure m shows a sketch of the embedding 
scheme as described above. During the self-consistency 
cycle we compute the particle number N^, a quantity 
that is obtained from the embedded Green’s function via 

= -^Tr J (28) 

To ensure particle number conservation, we need to up¬ 
date the electron chemical potential every time we re¬ 
ceive a converged self-energy from the inner-loop. For 
the present test cases, we found that the change in the 
chemical potential is relatively small, as demonstrated in 
Fig. [3] for bulk silicon (Si). However, we expect it to be 
more important for metallic systems. 



FIG. 3. Typical change in the chemical potential during the 
self-consistency cycle for a bulk Si calculation. Convergence 
is reached after 5 iterations of the outer loop. Already at 
the second iteration the chemical potential is close to its con¬ 
verged value. 


E. Total energy calculation 


Once self-consistency in the embedding scheme is 
reached, we can compute the total energy of the entire 
systems (the embedded cluster plus the environment) us¬ 
ing the converged lattice Green function Gy (k, w). The 
actual total-energy expression depends on the chosen 
methodology used in the embedded region. For hybrid 



FIG. 4. The embedded (Eq. (USD and the on-site (Eq. ([7|)) 
Green’s functions define two Dyson equations that form two 
nested loops. The two loops are iterated until self-consistency 
is reached. 


density functionals we have 

^^tj,(k)ny(k) 
k i,j 

+ E,,[n]+E^^[n.,], (29) 

where 

(30) 

R 

is the matrix form of the kinetic energy operator, E^s 
the electrostatic (Hartree plus external) energy, and E^l^ 
the XG energy. In Eq. (l29l) is the k-dependent global 
density matrix, 

«y(k) = y'^G(f(k,*u;)e—(31) 

and n is the electron density obtained from ny (k). We 
note that Eq. (l29l) is the exact total-energy expression for 
the hybrid density functional, and the only approxima¬ 
tion is that the density matrix (k) (and hence electron 
density n) is obtained from the RDMFE scheme and not 
from a periodic hybrid functional calculation. 

However, Eq. (l2^ cannot be directly applied since 
evaluating E^"^ as a functional of the k-dependent den¬ 
sity matrix (k) requires the computation of the exact- 
exchange energy for the entire periodic system, which 
is exactly what we are trying to avoid here. Instead of 
evaluating E^^[nij{k)\ in full, we thus only compute the 
change of E^^ with respect to the local or semi-local 
(EDA or GGA) energy in the embedded region. This is 
the main approximation of our approach, which is con¬ 
sistent with the spirit of the local self-energy correction 
in the RDMEE scheme, and is suggested by the near¬ 
sightedness of the XG energy of a bulk system (although 
the exact-exchange energy is probably not the most near¬ 
sighted self-energy we could have chosen) 


phyb __ 
"^tot 


Ni 


BZ 
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The Hartree energy depends on the electron density 
in a highly non-local way and it is questionable if a lo¬ 
cal treatment can be applied to the Hartree energy at 
all. Therefore, for simplicity, we omit possible changes in 
the Hartree and the external energy for now, assuming 
that the electron density given by the local or semi-local 
approximation is already sufficient. 

Finally, we are left with the kinetic energy term which 
also changes when irroving from local or semi-local to hy¬ 
brid functionals. For consistency, kinetic and XC energy 
should be taken together. In our scheme, we evaluate the 
changes of the kinetic and XC energy caused by the local 
self-energy correction within the embedded region. 

Based on the above considerations, we propose the fol¬ 
lowing approximate total-energy expression for embed¬ 
ded hybrid functional calculations 


phyb ^ rpKS 
^tot ~^tot 





1 pKSr KSi 

J ~ ^XC \ 


(32) 


where n™* is the embedded density matrix as defined 
in Eq. m, and is the “on-site” density matrix of 
KS-LDA/GGA calculations. n^®(r) is obtained from the 
on-site KS density matrix 


nKS(r) = ^ , (33) 


and is thus restricted to the embedded region. 

For GW we can proceed in an analogous fashion 


j^GW _ 


N, 





+ Aes[n] 

k i,j 

1 pGW r^lati 

+ Axe [Eij \ , 

(34) 

where 




j^GW r/^lati 

^XC i^ij J 

1 ^ 1 


-iui0+ 

2 “ W.BZ 

k 

(35) 


following directly from the Galitskii-Migdal (GM) 
formulai^ Similar to the hybrid functional case, we will 
not take the full k dependence in [G)^*(k, iw)] into 
account. Instead we adopt the same philosophy as before 
and make a local approximation 


T7’' 


GW 




KS 




- nW) 


tot ' / 

^3 

+ Axe l^ij (*^)J - Axeffi J , 


KS r^KSi 


(36) 


where 


To summarise this part, in RDMFE the total energy of 
the entire system can in principle be obtained from the 
lattice Green’s function. However, in practice, approxi¬ 
mations are needed to make the problem tractable. The 
expressions for hybrid functional and GW calculations 
proposed above are consistent with the local nature of 
the self-energy approximation in RDMFE, but their per¬ 
formance needs to be checked in practical calculations. 
Future work needs to revisit total energy calculations in 
RDMFE. 


III. COMPUTATIONAL DETAILS 

We used tight FHI-aims integration grids for all our 
RDMFE calculations. Eor the embedded PBEh and HSE 
self-energy we used the tier 1 basis set. Eigure[5]shows 
the embedded PBEh total and cohesive energies with 
increasing basis set size. The lattice Green’s functions 
were represented on a logarithmic frequency grid with 40 
points. The total energy calculations for 2 and 8 atom 
silicon unit cells and the density of states (DOS) calcu¬ 
lations for 8 and 16 atom unit cells were performed on a 
4x4x4 k-point grid, which we increased to 10 x 10 x 10 
for DOS calculations in the 2 atom unit cell. GW cal¬ 
culations were performed in a tier 3 basis set with 40 
frequency/time points in the inner loop and the same 
number of k-points as for PBEh. The linear mixing pa¬ 
rameter A was fixed to 0.5, which gave reasonably fast 
convergence. The periodic PBE and PBEO reference cal¬ 
culations were performed using the tier 1 basis set and a 
12 X 12 X 12 k-mesh. 

Densities of states are obtained in two different ways. 
In PBEh and HSE we obtain a self-energy that de¬ 
fines a converged k-dependent embedded Hamiltonian 
via the converged lattice Green’s function once the self- 
consistency cycle is converged. Eor the PBEh self-energy 
we can directly diagonalize the embedded Hamiltonian 
= H°(k) + ^t each k-point, which 

yields k-dependent eigenvalues and eigenstates. The re¬ 
sulting density of states (DOS) is n(e(k)) = X)'^(^(k) — 

ey(k)), where v labels the eigenstates of i7®™’’’®'^(k). To 
make the comparison to experiment easier, we introduce 
a Gaussian broadening 


exp 

to obtain the DOS iV(e) = ^ gi,(e). In this work we use 

a Gaussian broadening of cr = 0.2 eV. 

Eor the GW self-energy, the spectrum at each k-point 
is directly given by the Green’s function as 

A(k, uj) = -iTr{ImG^^*(k, w)}. (39) 

TT 


Ei^iGTt] = 11 ^Si°“(*a.)G^“‘’(*a.)e—. (37) 


To determine G*’^*(k,a;) on the real-frequency axis, we 
analytically continue the self-energy from the imaginary 











FIG. 5. Embedded PBEh total energy (upper panel) and 
cohesive energy (lower panel) for bulk Si with a 2 atom unit 
cell as function of the basis size. The energy zero is set at the 
value of the tier 1 basis. 


to the real axis. In practice, we fit a two-pole model, that 
has proven to work very well for the systems we tested, 
to each matrix element of the self-energ y 

where and /3„ are complex fitting parameters. We 
then evaluate eq. dMl) for real frequencies and solve 
Dyson’s equation for G'®'*(k,w). The spectral func¬ 
tion subsequently follows from a k-summation A{uj) = 
^ kl(k, w), which we convolute with Gaussians as 

k 

A{u}) = J duje A(uj), (41) 

with a broadening that we choose to be cr = 0.01 eV to 
obtain a DOS A(u)) that we can compare with experi¬ 
ment. 


IV. RESULTS 

Having introduced the concept of RDMFE and our 
implementation in the previous sections, we now turn 


to benchmark calculations for hybrid functionals, for 
which we have an independent, periodic reference in FHI- 
aims^. Then we present self-consistent GW calculations 
for which such a periodic reference does not yet exist in 
FHI-aims. We choose bulk Si as test system since it is a 
reliable and well studied reference case. 


A. Density of states and band structures 

We begin our benchmark tests by calculating the DOS 
at each iteration to investigate its evolution with each 
embedding cycle. Figure [5] shows the DOS at different 
iterations of the outer loop for a 2 atom unit cell of sil¬ 
icon. We observe that the largest change occurs at the 
first iteration when moving from PBE to our embedded 
PBEO DOS. For subsequent iterations the DOS changes 
are much smaller. 



FIG. 6. DOS comparison at each iteration of the main-loop. 
Gonvergence was achieved after 5 main-loop iterations. 


When comparing the converged, embedded DOS for 
the 2 atom unit cell with the periodic PBE and PBEO 
DOS shown in Fig. [71 we observe that the band width 
and band gap are larger than in PBE and are closer to 
the PBEO reference. When increasing the unit cell size 
to 8 and 16 atoms (see Fig. the difference between 
the embedded DOS and the periodic PBEO DOS reduces 
systematically. The resulting RDMFE band gaps for dif¬ 
ferent unit cell sizes are compared with the PBE and the 
PBEO values in Tab. HI With increasing unit cell size, the 
band gap increases and approaches the PBEO value. 

Next, we consider the band structure for the 2 atom 
unit cell shown in the upper panel of Fig. [HI We see the 
same trend as for the DOS: the band gap and the band 
width approach PBEO and so do the bands in general. 
However, at some high symmetry points the degeneracy 
of certain bands is lifted. The origin of this degeneracy 
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PBE PBEO 

RDMFE@PBE0 

2 atoms 8 atoms 16 atoms 

experiment (at 300K)^ 

band gap [eV] 0.68 1.85 

1.2 1.257 1.569 

1.12 


TABLE I. Comparison between PBE, PBEO and RDMFE for different unit cell sizes for the indirect band gap of silicon. The 
experimental valued is shown for reference. 



FIG. 7. Comparison of the embedded DOS, the periodic PBE 
and the periodic PBEO DOS for a 2 atom nnit cell. 




FIG. 8. Gomparison of the embedded DOS, with the periodic 
PBEO DOS for an 8 atom unit cell (npper panel) and a 16 
atom unit cell (lower panel). 


lifting is the break of the crystal symmetry that we in¬ 
troduce with the local self-energy. It is a well known 
artefact and has been discussed extensively in the con¬ 
text of cellular and cluster DMFT— The local self¬ 
energy simply does not “know” about the symmetry of 
the crystal and can therefore not enforce it. The solution 
to the problem is then obvious: the approximation of the 
locality of the self-energy needs to be improved. If the 
self-energy would extend over a larger region (i.e. super¬ 
cell) it would acquire more information about the crystal 
symmetry. Then the degeneracy splitting should reduce. 
In the two lower panels of Fig. we present an unfolded 
band structure^^ for the 16 and 32 atom unit cells. We 
indeed observe a reduction in the splitting for both the 
16 and 32 atom unit cells. However, while the degener¬ 
acy is fully restored for some high symmetry points, it is 
still broken for others such as the X and Z points. 

We will now turn to the GW spectra. The total spec¬ 
tral function for bulk Si with 2 atoms in the unit cell at 
the 1st iteration is shown in Fig. (TUI Since Dyson’s equa¬ 
tion has been solved once, this spectrum is not equivalent 
to perturbative GqWq spectra and we would expect to see 
plasmon satellites. The spectrum shows a broad peak be¬ 
tween -40 and -30 eV, which has been identified as plas¬ 
mon satellite^Si^ Such satellites are completely absent 
in KS band structures or in GqWq, because GqWq only 
corrects the KS states and does not yield new states. 
The energy range of the RDMFE satellite agrees well 
with previous periodic GW calculation s ^!^ and demon¬ 
strates that our dynamic, local RDMFE framework can 
capture non-local phenomena such as plasmon satellites. 
For sc GW the converged DOS is also shown in Fig. [TUI 
As demonstrated by Holm and von Barth^ for the elec¬ 
tron gas, full self-consistency in G and W leads to a dete¬ 
rioration in the GW spectral function due to the neglect 
of vertex corrections. Thus, the fact, that the plasmon 
satellite disappears at self-consistency is not surprising. 

We obtain a band gap of ^0.9 eV for the two atom 
unit cell, which is close to the experimental value of ^1.12 
eV— . This comparison together with the one between the 
indirect band gap from our calculation and experiment!^ 
are presented in Tab. im 


B. Embedded total energies: 

For scGW we currently do not have a periodic refer¬ 
ence to compare to, as alluded to before. We can, how¬ 
ever, construct another test case and benchmark against 
our scGW implementation for finite systems^i^, where 
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band gap 

RDMFE@scGW periodic sc GW-^ 

QPscGW^ experiment (at 300K) 

direct (Tisc) [eV] 
indirect (Eg) [eV] 

3.7 — 

0.9 1.55 

3.47 3.4^^ 

1.25 I.I 3 S 


TABLE 11. Direct and indirect band gaps as calculated from RDMFE for the GW self-energy. Comparison is made with the 
periodic scGW work of Kutepov et al— and the quasi particle self-consistent GW calculation of Kotani et al^ and experiment. 


the total energy was computed from the Galitskii-Migdal 
formula^i. We achieve this by considering the molecular 
limit of a unit cell, i.e., the limit of an isolated unit cell 
with a lattice constant of ~20 A. The benchmark results 
for He, H 2 and Na 2 are presented in Tab. im which shows 
the XC components that enter the total energy as given 
by Eq. (1571) . molecular scGW XC self¬ 
energy and the corresponding Green’s func¬ 
tion at convergence. is local XC self-energy and 

the embedded Green’s function both obtained 
at convergence of the RDMFE cycle. Table HlD illustrates 
that the components entering the embedded total energy 
agree almost to the meV level with the corresponding 
components from the finite systems sc GW calculation, 
which demonstrates the reliability and robustness of our 
implementation. 

We then investigated the convergence of the total en¬ 
ergy with respect to the increase of the unit cell size 
for RDMFE PBEO and sc GW. For embedded PBEO, we 
performed calculations for bulk Si up to 32 atoms in 
the unit cell, whereas for GW we considered bulk He 
in the fee structure up to 64 atom unit cells. To reach 
larger systems, a full parallelization of our implementa¬ 
tion would be required. The upper panel of Fig.fTTlshows 
the comparison of our embedded PBEO cohesive energy 
with the periodic PBE and the periodic PBEO energy. 
We also include a third reference in which we added the 
kinetic and XC energy of a PBEO calculation to the PBE 
energy, which most closely resembles our RDMFE ap¬ 
proximation. We see that with increasing unit cell size 
the embedded cohesive energy approaches the periodic 
PBEO value, but then dips below. This is not surprising 
since our embedded cohesive energy does not account for 
changes in the electrostatic energy. Instead, the RDMFE 
curve approaches the PBEO reference value from which 
the electrostatic change has been removed. However, the 
convergence to the periodic limit is relatively slow. This 
can be related to the long range nature of the HE exact- 
exchange as we will show later on using a range separated 
self-energy (see the discussion of Fig. [HD. For the GW 
self-energy, however, the total energy seems to converge 
much faster and only changes in a very small range. This 
is shown in the lower panel of Fig. Bll 

To better visualise the interplay between locality and 
unit cell size in our scheme we also present results of 
the BSE range-separated exact-exchange self-energy. We 
vary the range separation parameter to model different 
degrees of locality, but keep the percentage a of exact 
exchanged fixed. We then translate the range separation 
parameter into a radius i?sphere in real-space using the 


relation 7 = and determine the number of atoms 

that fit inside. We have considered range-separation pa¬ 
rameters 7 that correspond to spheres enclosing 2, 4, 8 
and 16-atom unit cells. For each 7 , the resulting embed¬ 
ded total energy is plotted in Fig. [T7]as a function of the 
size of the unit cell (upper panel). The lower panel shows 
the volume of the surrounding sphere for the different 7 
parameters and for the different unit cell sizes. We in¬ 
deed observe that the total energy converges faster with 
unit cell size, the shorter the range of the non-locality 
in the BSE self-energy. This proofs that RDMFE be¬ 
comes a viable option for self-energies, whose range only 
encompasses a few nearest atoms. In that sense, PBEO 
had been the toughest test, because its range is infinite. 

Finally, we briefly address cohesive properties. In 
Fig. [T5]we show the total energy of bulk Si as a function 
of the lattice constant for RDMFE@PBE0, scGW and for 
periodic PBEO. For RDMFE@PBE0 (upper panel) our 
calculations for a 2 atom unit cell already give a min¬ 
imum that is below the PBE one and very close to the 
experimental value of S.dSA^i (see Tab. IIVI) . However, for 
an 8 atom unit cell the lattice constant reduces slightly. 
For RDMFE@GW, the minimum for the 2 atom unit cell 
(lower panel) also lies below the PBE value and already 
agrees fortuitously well with the experimental value. Our 
RDMFE@GW lattice constant is slightly larger than 
that reported by a recent periodic self-consistent GW 
calculation^. We also performed a Birch-Murnaghan^ 
fit of the total energy curves to extract the bulk moduli 
for RDMFE@PBE0 and RDMFEOGW. The resulting 
values are reported in Tab. El 


V. DISCUSSION 

We have presented an embedding scheme for periodic 
systems that builds on DMFT. In our approach, the elec¬ 
tron interacting across periodically repeated unit cells is 
mapped onto an on-site problem, in which the electrons 
only interact directly in one unit cell, but are dynamically 
coupled to a periodic bath of electrons. The coupling be¬ 
tween the embedded system and the surrounding is con¬ 
structed naturally by means of Green’s functions. Due to 
its dynamic nature the bath can exchange electrons with 
the embedded region. Our embedding scheme is most 
suitable for systems with periodic boundary conditions, 
as the translational symmetry is automatically preserved. 
Furthermore, we transfer the non-locality of a self-energy 
into a frequency dependence, a concept that has previ¬ 
ously been explored by Gatti et alM- and in the spectral 
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Term 


He Hz Naz 


^GWj^embj ^ 1 J -27.173906 -17.503642 -760.630196 

EGw^QaaGWi^ = i / -27.181325 -17.505250 -760.632183 

^.GWjfjembj ^ 1 J^j.ioc(^^)gemb(^^)g-i^o+ -1.744578 -2.3053908 -2.457104 

^GW[(;jscGvr] ^ 1 J ^j.«:Gw^(.^)^<,cGvr(^^)g-iu;0+ -1.737910 -2.304200 -2.455707 


TABLE III. Exchange (X) and correlation (C) components of the RDMFE@scGIE total energy, as given by Eq. 

(EZJ in the limit of an isolated unit cell (large lattice constants): a benchmark against the standard sc GIT calculation for finite 
systems^, labeled by [G“°®'^]. All energies are in eV. 



RDMFE@PBE0 

periodic PBEO 

RDMFEOscGIF 

periodic sc GIF— 

experiment 

unit cell size 

2 atoms 

8 atoms 

2 atoms 

2 atoms 

— 

— 

lattice constant [A] 

5.45 

5.4 

5.43 

5.43 

5.39 

5.43^ 

bulk modulus Bq [GPa] 

95.14 

129.07 

84 

80.37 

100.7 

99^ 


TABLE IV. Bulk Si RDMFE equilibrium lattice constant and bulk moduli Bq for the embedded PBEO and GIF self¬ 
energies. Comparison is made with periodic PBEO performed with FHI-aims, the periodic sc GIF work of Kutepov et al^ and 
experiment 


density-functional theory of Kotliar et al^ We note that 
the only approximation introduced in our scheme is that 
the non-local XC coupling between neighboring unit cells 
(or computational supercells) is included only at the KS 
GGA level, and neglected in the more advanced (here 
hybrid functional or GW) treatment. In other words, 
the self-energy correction to the GGA XG potential is 
k-independent, an intrinsic feature of DMFT. 

We now compare our scheme to other embedding 
schemes. For the hybrid QM:MM approach a clear sepa¬ 
ration between the embedded region and the surrounding 
and the treatment of the boundary atoms is not always 
obvioua^^iSE. For systems, in which classical electrostat¬ 
ics dominate such as ionic or molecular solids, the sepa¬ 
ration between ions and molecules is natural. However, 
for covalently bonded systems it becomes more difficult 
to define the QM:MM partitioning. Thus, typically co¬ 
valent bonds are cut at the QM:MM boundary, which 
produces dangling bonds that need to be saturated. A 
multitude of models with different levels of accuracies 
have been developed to tackle these issues. One example 
is the ChemShell frameworb^SiSL that supports Hartree- 
Fock and hybrid functionals in the embedded region and 
that has recently been coupled to FHI-aims^. 

Another popular approach is ”Our own A-layer inte¬ 
grated molecular orbital molecular mechanics” (ONIOM) 
by Morokuma and coworkersj ^^'^^ ONIOM is a so-called 
extrapolative (or subtractive) scheme in which the total 
energy of the whole system is given by 

i^ONiOM = Arl — Aml + Amh, (42) 

where the RL refers to the real (or full) system at the 
lower level, ML refers to the model (or embedded) sys¬ 
tem at the lower level and MH labels the model system 


for the higher level theory. In constrat to the additive 
QM:MM scheme, ONIOM does not need an additional 
coupling Hamiltonian to describe the QM/MM intera- 
tion. When a QM/MM boundary cuts through a cova¬ 
lent bond, link atoms (mostly hydrogen atoms) are added 
to cap the unsaturated QM boundary for the model cal¬ 
culations. Even if it is common to use MM methods 
for describing the surroundings, the ONIOM scheme was 
recently extended to deal with two-layer two-QM em¬ 
bedding, ONIOM(QMl:QM2), where HF was used for 
the surroudings and the embedded model region is de¬ 
scribed by MP2 or B3LYP. The QM1/QM2 interactions, 
including electrostatic interaction, mutual polarization, 
and charge transfer, are described at the lower QM level. 

Our RDMFE scheme is distinctly different from the 
ONIOM(QMl:QM2) approach. First, the RDMFE 
scheme is formulated in terms of Green’s functions, 
whereas ONIOM is based on a partition of total ener¬ 
gies. As such, spectral properties come out naturally 
from RDMFE, while the evaluation of total energies is 
more involved, as discussed in Sec. Ill El The opposite is 
true for the ONIOM(QMl:QM2) scheme. Second, within 
RDMFE, the effect of the environment is encoded in the 
bath Green’s function that describes an electron reservoir 
with which the embedded cluster can exchange electrons 
freely. In other words, the electronic states in the embed¬ 
ded system are not forced to localize within the cluster, 
but are allowed to delocalize into the surrounding sys¬ 
tem. Thus, dangling bonds pose no conceptual problem 
and boundary effects are not significant since they di¬ 
minish quickly as the size of the cluster increases. In 
contrast, in ONIOM(QMl:QM2), like in most other em¬ 
bedding schemes in computational chemistry, link atoms 
are needed to saturate the dangling bonds when chemi- 
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Bulk Si with 2 atom unit cell 


embedded band structure - periodic PBEO 



Bulk Si with 16 atom unit cell 



Bulk Si with 32 atom unit cell 



FIG. 9. Upper panel: the embedded band structure for bulk 
Si compared to the periodic PBEO one. The local self-energy 
breaks the translation symmetry and the degeneracy gets 
shifted at some hight symmetry k-points. Middle panel: the 
embedded unfolded band structure for the 16 atom unit cell. 
The degeneracy shifting gets reduced compared to the 2 atoms 
case. Lower panel: the embedded unfolded band structure for 
the 32 atom unit cell. Also here the degeneracy is restored in 
most of the high symmetry points. 



FIG. 10. Gaussian broadened (with broadening cr = 0.01 eV) 
quasiparticle spectrum for the GW self-energy at 1st itera¬ 
tion (red curve) and at self-consistency (blue curve). Only 
occupied states are shown. 


cal bonds are broken. Therefore ONIOM(QMl:QM2) is 
most appropriate for describing systems with localized 
electrons, whereas RDMFE has no problem in dealing 
with delocalized electrons, especially metallic systems. 
Third, RDMFE, as is formulated right now, is only appli¬ 
cable to periodic systems that are relevant to solid state 
physics, while ONIOM(QMl:QM2) is most suitable for 
describing molecules and clusters that are of interest to 
chemical and biological applications. 

Addressing the problem of CO adsorption on 
Cu(lll)i^, Hu, Reuter and Scheffler— developed a 
cluster extrapolation scheme that is based on perform¬ 
ing a cheap (LDA/GGA) calculation for the periodic 
system then correcting the resulting total energy by 
AExc = [LDA/GGA] - “better”], where 

^duster[ ldA/GGA] and “better”] are the clus¬ 

ter XG-energy parts of a cluster calculation with the 
cheaper (LDA/GGA) and the ’’better” theory respec¬ 
tively while the cluster itself has been cut out from the 
periodic system. Increasing the cluster size, they could 
then show that the correction AExc converges for rel¬ 
atively small cluster sizes (^16 atoms) and thus much 
faster than [” tetter” ] alone. This cluster extrap¬ 

olation concept is similar to that of ONIOM(QMl;QM2) 
described above, but link atoms were not used for the 
cluster calculations. 

Whitten and coworkers also approached molecular ad¬ 
sorbates on metal surfacesi^ They developed an embed¬ 
ding scheme that builds on identifying a localized sub¬ 
space that has maximal exchange overlap with the va¬ 
lence orbitals of the atoms within and bordering the ad¬ 
sorbate. The localized subspace is then solved using the 
Cl method, for a fixed Coulomb and exchange potential 
constructed from the localized orbitals. However, the 
approach mimics the real periodic system using a large 
cluster of atoms, which fails in describing the system ac- 
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FIG. 11. Upper panel: the embedded PBEO cohesive energy 
for bulk Si for increasing unit cell size. Lower panel: the 
embedded GW total energy for bulk He for increasing unit 
cell size, where the blue region indicates a change in the range 
of ~ 20 meV. Both curves are referenced to the PBE total 
energy. 


curately. Additionally, no systematic cluster extrapola¬ 
tion has been studied. In a similar spirit, Huang and 
Carter developed a density-functional-based embedding 
scheme The scheme relies on the fact that the density 
is additive, i.e., that the total electron density can be 
partitioned into the density of the embedded region and 
the density of the embedding surrounding. Proceeding 
as such, allows the definition of an embedding density- 
potential, that is a functional of the total and the em¬ 
bedded density. Adding this potential to the embedded 
Hamiltonian and solving the resulting KS Schrddinger 
equation self-consistently leads to the desired embedded 
density. For the embedded region, correlated wave func¬ 
tion methods are typically used, while for the embedding 
potential the optimized effective potential method or ki¬ 
netic energy density functionals are employed. Moreover, 
due to the static nature of the embedding potential no 
dynamical methods can be used to describe the embed- 
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FIG. 12. Upper panel: the embedded HSE total energy for 
bulk Si for increasing unit cell size and different screening pa¬ 
rameter 7 . The range separated total energy converges faster 
than the PBEO one (black curve). All the embedded HSE 
total energies are shifted to the 2 atom unit cell value of the 
embedded PBEO curve to allow better comparison. Lower 
panel: The change of the sphere radius and screening param¬ 
eter 7 (on a reciprocal scale) with the volume of the sphere 
surrounding the 2, 4, 8 and 16 atom unit cell. 


ded region, which limits the applicability to ground state 
properties. 

For point defects in semiconductors, Schefher et al^ 
devised a self-consistent Green’s function method to com¬ 
pute the change in density induced by the presence of the 
defect. They considered this change as being a perturba¬ 
tion to the perfect crystal and solved the resulting Dyson 
equation self-consistently. Using the fact that defects are 
well localized in real space, they correct the Hellmann- 
Feynman force of the perfect crystal - calculated with 
force fields - by a contribution containing the change 
in density due to the defect - calculated with KS-DFA. 
They showed that the resulting Hellmann-Feynman force 
is comparable in accuracy to a full DFA calculation. 

Finally, it is also worth mentioning, that we perform 
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Lattice constant [A] 



Lattice constant [A] 

FIG. 13. The embedded total energy for bulk Si for different 
lattice constants: RDMFEQPBEO is shown in the upper and 
RDMFE@GIT in the lower panel. 


fully self-consistent GW calculations^^^ in our scheme, 


which is conceptually different from the so called quasi¬ 
particle self-consistent GW concep t^^'^^ (QPscGtF). In 
QPscGIF a series of GolPo calculations is performed. In 
each iteration the “best” Go is determined that most 
closely resembles the GW Green’s function of the cur¬ 
rent cycle. In practice a static, non-local potential is con¬ 
structed that approximates the GqWq self-energy. This 
non-local potential defines a new non-interaction Hamil¬ 
tonian Hq that produces a new input Green’s function 
Gq. Since the QPscGIF concept also requires the calcu¬ 
lation of the full non-local GqWq self-energy our expec¬ 
tation is that it will be easier to go beyond GW in our 
RDMFE framework. 


VI. CONCLUSIONS 

We have presented an embedding scheme for peri¬ 
odic systems based on Green’s functions in the DMFT 
framework, which maps an infinite periodic system to a 
single-site (single unit cell) problem coupled to an elec¬ 
tronic bath that needs to be determined self-consistently. 
Our RDMFE Green’s function mapping allows a nat¬ 
ural definition of the embedded region and defines a 
self-consistency loop which, at convergence, yields self- 
consistent Green’s functions. The coupling to the sur¬ 
rounding is of dynamical nature enabling electron ex¬ 
change between the embedded region and the surround¬ 
ing. We showed that our scheme produces densities of 
states and total energies that converge well with increas¬ 
ing size of the embedded region. We also demonstrated, 
that the main features of the “better” theory are rapidly 
captured within our scheme; for example the plasmon 
satellite already appears in RDMFE@GIF calculations 
for 2 atoms in the Si unit cell. RDMFE is therefore a 
promising embedding scheme, that has the potential to 
make sophisticated and computationally expensive first- 
principles theories available for periodic systems. 
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